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Abstract. The implementation of global non-analytical fit is exemplified by its 
application to the spectral analysis of the complex lineshape of P^S^CaC^Os+s 
at (7r,0). It deals with properly a multi-level globality in the fitting parameters 
and efficiently the non-analytical evaluation of the fitting function, thus exhibiting 
a potential applicability to a wide range of systematic analysis task. 

PACS numbers: 74.72.Hs, 71.10.Ay, 79.60.Bm, 02.60.Ed, 75.60.-d 
Angle-resolved photoemission spectroscopy (ARPES) has greatly advanced our 
microscopic understanding of the high-temperature superconductivity by presenting 
high-precision measurements of various dependence P^. Among them, for example, the 
doping and temperature (T) dependent measurement on Bi 2 Sr 2 CaCu208+5 (Bi2212) 
revealed the existence of a coherent component in the lowest-lying binding energy 
position which is closely related to the superconducting transition [2]; the momentum 
and photo energy {hv) dependent investigations identified the long-sought bilayer band 
splitting effect in the very overdoped (OD) Bi2212|Sj and argued for its persistence 
across a large doping range of the sample|lj. Due to the spectroscopic nature of ARPES, 
the intrinsic single-particle spectral function is masked by the matrix element effect jS] 
from a direct detection of the multi-component lineshape, for example, the bilayer- 
split peak- dip-hump lineshape of the OD65 Bi2212 (with a T c =65 K). In Ref.jH], a 
combined analysis of the hv and T dependence was achieved for a reliable extraction 
of the intrinsic spectral information. To naturally account for the physical constraints 
imposed by the intimate connection of the two dependence of data, a global fitting 
program was developed, to our knowledge, for the first time in the spectral analysis of 
complex lineshape with simultaneous consideration on various dependence. 

In this paper, we revisit this systematic spectral analysis on the 35 energy 
distribution curves (EDC's) [Fig. l(a)-(c)] by detailing the implementation of the global 
fitting procedure. It deals with properly a multi-level globality (defined below) in the 
fitting parameters and efficiently the non-analytical evaluation of the fitting function, 
thus exhibiting a potential applicability to a wide range of systematic analysis task. 

f Electronic mail: ruihuahe@stanford.edu; Present address: Department of Applied Physics, Stanford 
University, Stanford, California 94305 
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Figure 1. (color) Results of the global non-analytical fit on the T-dcpcndcnt (a) (sec 
also [d] for an enlargement of the T=30 K curves) and /^-dependent ARPES spectra 
at (it, 0) of (Pb)-OD65 in the superconducting state at T=10 K (b) and normal state 
at T=80 K (c). The AB hump lies lower in energy than the BB hump. The width 
difference of the peak is due to resolution variations from 10 to 18 meV at different 
hu's. 

For clarity, we recap the fitting function used for the lineshape modelling in Ref. 
jH] with only a slight modification described below. 

4 

I(u, T, hv) = J (T, hv) ■ [QT J a (u, T, hu) ■ f{u, T)) ® R(u, T'(hu)) + B(u, T)\ + h(T, hu), 

a 

, where / is the Fermi function, R the /izz-dependent resolution Gaussian, B the In- 
dependent empirical background function[H], Iq and I\ the EDC-specific linear intensity 
coefficients. The summation is over the spectral intensity of anti-bonding band (AB) 
hump, bonding band (BB) hump, AB peak and BB peak, given by, respectively, 

J a h{v, T, hv) = M ah (hv) ■ C a {T) ■ A h (u, T, a, e ah ), 
J bh (uj, T, hv) = M bh (hv) ■ C b {T) ■ A h (u, T, a, e bh ), 
J ap {u, T, hu) = M ap {hv) ■ (1 - C a (T)) ■ A p {u, r a (T),e ap ), 
J bp (cu, T, hu) = M bp {hv) ■ (1 - C b {T)) ■ A p {u, T b (T),e bp ) 
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, where M is the (squared) matrix element dependent only on hi>, C the remaining 
spectral weight in the hump with its T dependence to be fitted, Ah and A p the self- 
normalized spectral functions for the hump and the peak, respectively, given by 

A ( T \ £IEV,T)| 



u - Jco 2 + A SC (T) 2 ) 2 + E>,T) 2 ' 



We modified the empirical form of the imaginary part of self-energy into £"(u;,T) 



(auj) 2 + ((3T) 2 + (\, where lifetime broadening associated with the impurity level of 
samples is reflected by ( and that with doping level by a and (3. They are assumed to be 
independent of hv and T. £ is a self-normalization factor given by Pfj^ A h {ui)duj = 1 [T] . 
T(T) is a T-dependent linewidth of the peak. For the bilayer-split bands in the 
superconducting state, a realistic SCG opening introduces an effective T-dependent 
energy shift from their T-independent renormalized band energy positions (e a h and £&/j) 
while the energy position of the peak (e ap or Sb p ) is subjected to the fit. 

In order to retrieve reliable quantitative information from a mult i- dimensional fit 
to a set of interconnected EDC's, many physical constraints are involved and strictly 
followed. For example, in the the /^-dependent set, all functions of T (C, for example) 
are shared variants to all EDC's at the same T while the functions of hu (the matrix 
elements, for example) are locally specified, and in the T-dependent set, vice versa. In 
other words, for a given fitting parameter which is subjected to a constrain, it is required 
to have the same value for a specific group of EDC's (forming a so-called EDC's subset 
for this parameter). Thus, the total number of its possible values across the whole EDC 
set is smaller than the set size, 35. We define this number the globality level (Globality 
for this fitting parameter. Particularly, the parameters with Globality = 1 are in the 
top globality level, for example, Eh, a and those kept fixed during the fit, such as the 
parameters for the BB peak component which is not considered in the three-component 
fit exemplified below; the lowest globality level (Globality = 35) is owned by those local 
parameters with their values specific to each EDC, such as I . 

Technically, our fitting function contains some non-analytical features, convolution 
by the resolution function and integration in the self-normalization of spectral function. 
The conventional point-by-point evaluation of the fitting function in the fitting iteration 
fails to guarantee an efficient optimization to such a large data set. We developed our 
fitting program in WaveMetrics Igor Pro 4.0[7J. After implanting the fast algorithm for 
non- analytical fitting |Hj in the two- level global fit procedure by Igor (where only the 
top and the lowest globality levels are concerned), we further establish its multi- level 
global fitting capability. We illustrate below the 5-step implementation of our global 
non-analytical fit based on the 3-component model (neglecting the Jb p term in the first 
equation) on the hv and T dependent lineshapes of (Pb-)OD65[6j. 

Step 1 : Include in the fit the 9 normal state (@80K) then the 9 superconducting state 
(@10K) /^-dependent EDC's on Pb-OD65 (in an ascending order sorted by the 
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Table 1. List of fitting parameters 



Para No. 


1 




2 


3 4 


Parameters 


h M ah 




M bh 


Map M bp 


Globality b 


35 10(2x9 + 17x1) 


10(2 x 9 + 17 x 1) 


10(2x9+17x1) 1 


Fixation 


X X 




X 


x arb. a 


Para No. 


5 


6 


7 8 


9 10 11 


Parameters 


C a 


c b 


£ah £bh 


£-ap &bp & 


Globality 


17(10 x 2 + 1 x 15) 


1 


1 1 17(10x2 + 1x15) 1 1 


Fixation 


X 


1 


X X 


x arb. x 


Para No. 


12 


13 


14 15 


16 17 


Parameters 


r a 


r b 


5E F d h 


C 


Globality 


2(18 x 1 + 17 x l) c 


i 


1 35 


1 2(18 x 1 + 17 x l) c 


Fixation 


X 


arb. 


x 


X X 



a arbitrary value; 

b the number without (or outside) the bracket: Globality (number of subsets); expression inside the bracket: size of 

each subset (number of EDC's contained) X number of subsets of the same size; 
c to assume that the quality difference between Pb-OD65 and OD65 used in the two dependent investigations should 

be reflected only in the quasi-particle lifetime; 
d a slight correction in energy for Ep of each EDC to account for the uncertainty in Ep alignment. 

hu values, respectively). Append to the EDC's set the 17 T-dependent EDC's on 
OD65 using He-I (in a descending order sorted by the T values). The EDC's are 
labelled from to 34 and concatenated sequentially in X and Y scale, respectively, 
to form a huge "EDC". 

Step 2 : Discriminate the fitting parameters of the top globality level from the otherwise 
(of intermediate globality or the lowest globality). Then fix those parameters if 
they are required to be invariable during the fit. We summarize in Table 1 the 
properties of all the fitting parameters involved. 

Step 3 : Build a global fit parameters set, which includes sequentially, the 9 parameters 
with Globality = 1, the remaining 9 parameters with Globality > 1 (9 x 35 = 315, 
totally) grouped in the number order of the EDC in the EDC concatenation to 
which they belong. Hence, we end up with another index to locate the parameter 
in a global manner, called GL Para No. hereafter, from which we can set up a 
bidirectional mapping to the original Para No. (in Table 1) for a specific EDC. 

Step 4 : The global fit problem on a set of EDC's has been converted into a fit on the 
EDC concatenation with the global fit parameters set. Having in mind the unique 
calling pattern of the fitting function in Igor curve fitting routine jH] and the time- 
consuming non-analytical evaluation in the fitting function (see Ref. jH]), a direct fit 
is practically unfeasible, which asks for a function calling time up to (M+l) x N (M 
and N denote the total GL Para No. and point numbers of the EDC concatenation, 
here 324 and ~ 5400, respectively) for a single iteration! The introduction of a fast 
algorithm greatly reduces the function calling time down to (M+ 1) x L (L denotes 
the total EDC number included in the concatenation, here 35) at the expense of 
a size increase in the RAM which is used for the storage of the trial fit curves |Hj. 
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Figure 2. Truth tables as black-and-white plots (black: True, white: False), where the 
true pertinence of elements in the global fit parameters set to each EDC is recovered for 
a) a two-level global fit or b) a multi-level global fit in the example problem. The arrows 
identified by the Para No. locate the corresponding features of the fitting parameters 
with intermediate globality levels. For example, the three arrows from "f" indicate 
M a h has two independent values for the hv and T dependence subset, respectively, one 
for the EDC pairs taken at the same hv and the other for all the EDC's taken with 
He-I (the so-called parameter degeneracy of M a h). A horizontal line segment with 
unit length stands for a parameter with Globality = 35 while a full-scale line segment 
for the one with Globality — I. Note that the patterns may vary dramatically with 
different orders the EDC and parameters of the sets would be arranged in. 



We can further cut down the time to (K + 1) x L (K denotes the total Para No., 
here 18) by exploiting the fact that in a global fit a large number of the global fit 
parameters are actually not related to the function evaluation for a given EDC, 
which needs to be reduced from the EDC concatenation during each iteration of 
the fit. This is because only K out of M parameters in the global fit parameters 
set comes from a single EDC. 

To pick out properly those relevant parameters which really contribute to the 
function update on a given EDC, a truth table helps, which is convenient to build at 
the end of Step 3. In Fig. E[a), we show this binary table as a black-and-white plot. 
The bottom black wall and black bricks in the diagonal clearly indicate the order 
in which the global fit parameters set is organized and, especially, the EDC-specific 
grouping nature of the Globality > 1 parameters following the Globality = 1 ones 
(reminiscent of Step 3). 

Step 5 : The global fit involves only two levels in Step 4. The task of this step is to 
achieve a multi- level globality defined in Table 1. This only requires our making 
modifications on the truth table accordingly. Fig. Efb) shows a truth table 
modified to achieve the wanted globality configuration. In contrast to Fig. |2Ja), the 
introduction of intermediate globality levels to the global fit parameters set results 
in the parameters degeneracies for certain EDC's as well as, in return, some "void" 
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parameters, which seemingly have no effect on any EDC in the concatenation. To 
prevent Igor curve fit routine from "a singular matrix error" report, we simply 
fix them to arbitrary values during the fit. This modification doesn't alter the 
conservation in the column sum of truth values in the table (i.e., the total Para 
No., K), and thus neither the iteration efficiency. 

In Fig. ^1 the fit curves collapse well onto the experimental EDC's, reminiscent 
of Figs. 2-3 in Ref. 6J. Pay attention to the great contrast in the overall quality (or 
details) of the fit between the ones shown in Fig. 2 of Ref. [§J and ours, though on data 
sets generally twice larger in size. This suggests the local fit scheme ITU] with manual 
adjustments on the global fitting parameters is unlikely to yield a reliable result in a 
truly global sense. In contrast, our fully-automatic fitting routine can yield a robust jTTj, 
efficient and physically-constrained global fit for a large data set, which is frequently 
encountered in the spectral analysis of ARPES. 

Note that the algorithm proposed here is not limited to the application of the 
ARPES lineshape analysis only, but of a great applicability to a wide range of problems 
where a multi-level global analysis or/and non-analytical nonlinear fit is involved. 
For example, the analysis of hysteresis curves in the magnetism research can be 
a good case in point for its potential application. The angle- dependent hysteresis 
measurement contains the information of magnetic anisotropy which is currently 
retrieved only qualitatively due to the lack of a global fit on the angle-dependent curves 
set. Furthermore, in most cases, the single-domain coherent rotation model in the 
textbook is not sufficient to include the necessary physics to reproduce the hysteresis 
curve in calculation. By using spin dynamics |1 2 j as the non-analytical evaluation kernel 
in the global fit, the roles of collaborative spin motion and domain formation which are 
requisite for the hysteresis can be investigated quantitatively. Generally, our algorithm 
provides a useful and universal guideline to serve the advanced systematic analysis in 
experimental physics. 
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